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!2t (57) Abstract: The invention concerns a tomographic image reconstruction computation involving the breakdown of two-dimen- 
S sional sets of projections of said image by a wavelet decomposition method, wherein the set of projections is separated into pictures 
fS expressing respectively the set and the details at a higher frequency of the projections. The reconstruction is performed separately 
q on the pictures before any combination is achieved. Some properties of said decomposition, in particular in the Fourier frequency 



domain, enable substantial reduction of the amount of computations. 
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(57) Abrege : Un calcul de reconstruction d'images en tomographic implique la decomposition de jeux bidimensionnels de pro- 
jections de cette image par une miSthode de decomposition en ondelette, dans laquelle le jeu de projections est separS en imagettes 
exprimant respectivement l'ensemble et les details a plus haute frequence des projections. La reconstruction s'effectue separement 
sur les imagettes avant qu'une combinaison ne soit faite. Certaines propri£t6s de cette decomposition, notamment dans le domaine 
frequentiel de Fourier, permettent de require sensiblement la quantity de calculs. 
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PROCEDE DE RECONSTRUCTION D * UNE IMAGE A PART I R D'UN JEU 
DE PROJECTIONS PAR APPLICATION D'UNE TRANSFORMEE EN 

ONDELETTE 



5 DESCRIPTION 

Le sujet de cette invention est un procede 
de reconstruction d'une image a partir de jeux de 
projections de cette image et en appliquant une 
transformation en ondelette . 
10 Les procedes de tomographie consistent a 

examiner un objet inanime ou un etre vivant par un 
reseau de detecteurs mobiles qui prend une serie de 
vues en tournant autour de lui . Ces vues sont des 
projections de la propriete permettant d'expr inter 
15 1' image (normalement une absorption d'un rayonnement 
traversant 1" objet ou une scintillation d'un corps 
radioactif ingere par l'objet), c'est-a-dire des sommes 
de la propriete le long de lignes traversant 1' objet et 
definies par la collimation des detecteurs. Chaque 
20 detecteur mesure une projection de 1' image a chaque 
vue. Quand un nombre suffisant de vues et de 
projections a ete pris, on procede a 1' inversion des 
resultats pour obtenir la valeur de la propriete a 
chaque point de 1' objet ; cette inversion est 
25 comparable a 1' inversion d'un systeme d' equations de 
grandes dimensions et peut s'effectuer ouvertement par 
des methodes algebriques, ou, plus frequemment, par des 
methodes analytiques par lesquelles on applique des 
operations numeriques successives aux projections sans 
30 inverser directement le systeme. Un grand nombre de 
methodes existent, parmi lesquelles on citera celle qui 
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est decrite dans le brevet francais 2 615 619 qui est 
le premier de cette equipe de recherche, et le brevet 
francais 2 810 141 plus recent qui presente quelques 
ressemblances avec le procede qui sera decrit ici. On 
5 peut aussi citer un article de Grass paru dans Phy. 
Med. Biol, vol.45, p. 329, fevrier 2000. Ces operations 
aboutissent a ce qu'on appelle la retroprojection des 
mesures, c'est-a-dire au calcul de la valeur de la 
propriety prise a chacun des points de la ligne de la 
10 projection. 

II peut etre avantageux de travailler avec 
des resultats exprimes dans le domaine frequentiel de 
Fourier, comme le prouve le deuxieme document. Des 
transformees numeriques d' autre nature ont aussi ete 
15 employees . 

On se reportera aux figures 6, 7 et 8 pour 
une description concrete et schematique d'un procede de 
prise de mesures tomographiques . Une source F de 
rayonnement et un systeme de detection 3 sont mobiles 
20 le long d'un cadre annulaire 2 a des positions 
diametralement opposees, et un rayonnement conique 
originaire de la source F atteint le systeme de 
detection 3 apres avoir traverse l'objet 1 a etudier. 
La piece essentielle du systeme de detection 3 est un 
25 reseau 4 bidimensionnel de detecteurs 5. Des 
projections R de 1 ' image tridimensionnelle de l'objet 1 
sont mesurees par ceux des detecteurs 5 qui sont inclus 
dans un pourtour 15 de 1 '"ombre" de l'objet 1. Un 
nombre important de vues de ce genre est pris a autant 
30 d' angles 9 differents d' orientation du reseau 4 des 
detecteurs 5. II est frequemment considere, 
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convent ionnellement pour les calculs, un reseau fictif 
4 1 de ditecteurs 5 1 sur un plan de detection Pdet qui 
est paralldle au reseau 4 rgel et passe par le centre O 
du cadre 2. Des coordonn^es p et q sont d6finies pour 
5 identifier les d6tecteurs 13 d'aprds leurs lignes et 
leurs colonnes. Des calculs de rearrangement, courants 
dans I'art, permettent de transposer les mesures d f un 
reseau 4 quelconque au reseau fictif 4' et d'appliquer 
les algorithmes de reconstruction a celui-ci. 
10 Le brevet frangais 2 615 619 sera rappel£ 

ici puisqu f il expliquait en detail quelles operations 
num^riques, notamment celles qui sont dites de filtrage 
et de r€troprojection, permettent d'obtenir 1 1 image de 
l'objet & partir de ses projections ; des methodes 
15 differentes existent cependant . 

Le probl£me auquel on s f attaque ici est la 
reduction du temps ou du volume des calculs pour 
inverser le systdme des mesures. 

On sait qu'il s ! agit d'une des limitations 
20 les plus graves dans les proced§s de tomographie, et de 
nombreux proc6d§s nouveaux ont 6te concpus pour le meme 
objectif que l f invention, dont celui du deuxieme brevet 
cit6 . 

L'id6e exposee ici est d' exploiter les 
25 propri6t6s particuliSres d'une transform^ num£rique, 
dite decomposition en ondelette, des projections pour 
prevoir des emplacements de projections nigligeables ou 
insignif iantes et de ne pas appliquer les calculs 
d f inversion & ces emplacements. Des perf ectionnements 
30 off rent encore une plus grande reduction des calculs. 
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Sous sa forme la plus generale, I 1 invention concerne un 
proced4 de reconstruction d'une image a partir de jeux 
de projections de cette image, comprenant 

successivement : 
5 - une serie de decompositions successives 

des jeux de projections en ondelette dormant des 
imagettes des jeux des projections, comprenant des 
images d'une approximation (AA) et des series 
successives (Dd, Dh, Dv) de details homologues de 
10 chaque jeu, 

- dans chacune des series et successivement 
pour des imagettes de details ayant un notnbre croissant 
de points, une recherche de portions insignif iantes 
jug^es pour depourvues de contenu, et une recherche de 

15 portions homothetiques des portions insignif iantes dans 
chacune des imagettes de details qui suivent dans la 
s£rie ; 

- des retroprojections des imagettes des 
jeux de protection en imagettes de 1 1 image k 

2 0 reconstruire, avec une omission des retroprojections 

pour toutes les portions insignif iantes et toutes les 
portions homothetiques. 

et une combinaison des imagettes k 
reconstruire par une inversion de decomposition donnant 

25 ladite image. 

L'ordre des etapes (surtout les 
retroprojections precedant la recomposition) est 
essentiel pour obtenir les avantages de 1' invention. 

L'art anterieur comprend un exemple de 

3 0 reconstruction d'une image par une decomposition en 

ondelette (US - 5 953 388 - A), qui cependant est 
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appliqu€e 1& pour reconstruire une portion seulement de 
l 1 image, en exploitant la propriSte de "localite" de la 
decomposition, qui n'est guere sensible aux autres 
portions de l 1 image, qu'on peut ainsi n^gliger dans les 
5 calculs. Le brevet US-5 841 890-A traite d'un sujet 
analogue et considere divers aspects de la 
decomposition en ondelette appliqu6e k la tomographic. 
La decomposition ne sert pas a acc£l€rer des calculs de 
retroprojection pour obtenir une image d' ensemble de 
10 fagon plus econome en calculs. 

On doit enfin mentionner le document US-6 
148 110-A qui decrit un proc6de a masques num€riques 
analogue a celui de l 1 invention, mais seulement pour 
comprimer un signal d 1 image sans calcul de 
15 retroprojection. 

L f invention sera maintenant d^crite plus 
concretement et completement en liaison aux figures 
suivantes : 

- la figure 1 illustre une decomposition 
20 d 1 image en ondelette, 

- les figures 2, 3 et 4 illustrent certains 
aspects de 1' invention, 

- la figure 5 illustre un organigramme 
resumant une realisation complete du proc6d6, 

25 - la figure 6 illustre un dispositif de 

principe du proc6d6, 

- et les figures 7 et 8 illustrent plus en 
detail le proc§d6 et certaines rotations employees. 

Nous commencerons par exposer la 
30 transformation d'un signal par une decomposition en 
ondelette. Plusieurs modeles d'ondelettes existent, qui 
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ont en commun d'etre comparables a un filtre passe-bas. 
Le signal est separe en deux portions dont l'une, 
associee aux basses frequences, peut §tre tenue pour 
une approximation du signal, alors que l'autre, 
5 associee aux hautes frequences, exprime plutot ses 
details. Une propriete des ondelettes est que les 
portions peuvent contenir chacun une moitie des points 
du signal si bien qu'il n'y a pas de perte 
d' information par cette decomposition. La decomposition 
10 peut etre faite dans le domaine direct d' expression du 
signal ou dans le domaine de Fourier. 

Dans le cas de projections d'un objet 
examine sur un reseau generalement bidimensionnel de 
detecteurs, les projections peuvent etre groupees en 
15 jeux bidimensionnels d'apres deux de leur coordonnees 
(en general p et q sur les axes d'un reseau de 
detecteurs). Toutefois, comme la figure 8 l'illustre, 
on considere le plus souvent des jeux de projections 
rearrangees sur un reseau fictif de detecteurs. Dans 
20 l'exemple de la figure 8, on regroupe des projections 
Rx originaires d'un certain nombre de positions Fx 
successives de la source F du rayonnement, de maniere 
que les projections aboutissant a une meme colonne (a q 
constant) de detecteurs 5' du reseau fictif 4' soient 
25 originaires d'une meme position Fx, et aussi que les 
plans de projections soient tous paralleles jusqu'aux 
positions Fx : on a alors transforme le probleme d'une 
geometrie conique du rayonnement a une geometrie en 
eventails paralleles dont la resolution et plus simple. 
30 De plus, le reseau fictif 4' passe par le centre 0 de 



FEUILLE DE REMPLACEMENT (REGLE 26) 



WO 2004/036288 



PCT/FR2003/050089 



10 



rotation et appartient done au plan de detection Pdet, 
ce qui facilite aussi les calculs. 

L 1 invention pourrait encore etre appliquee 
a des rencont ructions de coupes de l'objet 1 a l'aide 
d'un reseau monodimensionnel de detecteurs (tous a la 
meme coordonnee p) . On pourrait encore appliquer le 
principe expose ci-dessus du rearrangement des 
projections d'une geometrie en eventail plan a 
geometrie parallele. 

Le traitement des mesures rearrangees se 
fait en suivant les lignes de detecteurs 5 du reseau 
fictif 4', successivement pour tous les points. 

Un signal peut etre decompose 
successivement en ondelette pour donner plusieurs 
15 niveaux de resultats. Les nouvelles decompositions ne 
concernent que la portion de 1' ondelette qui donnait 
1' approximation du signal, la portion ou les portions 
qui donnaient les details etant preservees. 

Prenons comme exemple d'objet d'une 
decomposition en ondelette une image formee de cinq 
cercles dont un cercle exterieur et quatre cercles de 
diametres differents, tous inscrits dans le premier. 
L' image de depart comprenait n x n points, et si on 
applique deux fois une decomposition en ondelette de 
cette image selon le principe precedent, le resultat 
est donne a la figure 1. 

La decomposition en ondelettes de 1' image 
donne un ensemble d'imagettes, dont trois sont plus 
importantes que les autres, comprenant n/2 x n/2 points 
chacune, et correspondent aux details horizontaux, aux 
details diagonaux et aux details verticaux de 1' image 



20 



25 



30 
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de depart a grande echelle ; elles sont notees 
respectivement Dhl, Ddl et Dvl . Les details horizontaux 
de 1' image sont obtenus a partir des projections 
d' angles 6 proche de O ou %, les details verticaux a 
5 partir des projections d' angles 6 proche de n/ 2 ou 3rc/ 2 , 
et les details diagonaux a partir des projections 
d- angles interm6diaires avec les conventions de la 
figure 1. Le reste de 1' image est compose de quatre 
imagettes comprenant n/4 x n/4 points chacune et dont 
10 trois sont des imagettes de details horizontaux, 
diagonaux et verticaux a plus petite echelle, notees 
Dh2, Dd2 et Dv2 f alors que la derniere imagette est une 
approximation notee AA de 1' image de depart. Si on 
appelle <|> et \|f les fonctions de decomposition d'une 
15 image ou d'une imagette en ondelette, la fonction <|> 
dormant 1' approximation et la fonction \|f les details, 
les fonctions a appliquer a 1' image de depart pour 
obtenir la decomposition de la figure 1 sont donnees 
par le tableau I . 
20 AA : <|> (Xi) <|> (x 2 ) <J> (-j) <t> <y> 

Dv2 : <|> (x x ) <|> (x 2 ) * (f) V (y), 

Dd2 , <> * (x 2 ) ¥ (^) V (f >, 

Dh2 : <> (X X ) $ (x 2 ) 4> (^) 4> (f) 

DVl : <j> (Xi) ¥ (X 2 ) 

25 Ddl : 4* (Xi) ¥ (X 2 ) 

Dhl : ¥ (Xi) 0 (X 2 ) 

L' invention consiste a effectuer la 
retroprojection sur chacune des imagettes des jeux de 
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projections decomposees en ondelette et a combiner les 
imagettes ret roproje tees par une inversion de la 
disposition en ondelette pour obtenir 1' image 
recherch6e. La decomposition en ondelettes est propice 
5 a diverses simplifications qui accelerent beaucoup les 
calculs. Ces simplifications sont faites entre la 
decomposition et la combinaison. 

La premiere d' entre elles concerne des 
filiations qu' on peut etablir entre des details 
10 homologues a des echelles differentes. On considere 
pour cela des series d' imagettes dormant des details de 
mime nature. La figure 2 (qui^ illustre une 
decomposition d'une image ressemblant a celle de la 
figure 1, mais a trois niveaux de decomposition) 
15 illustre, pour les trois imagettes des details 
horizontaux Dhl, Dh2 et Dh3, des portions homologues 
Jl, J2 et J3 qui occupent la meme position et la meme 
superficie relative sur chacune de ces imagettes et se 
deduisent ainsi l'une de 1 ' autre par une homothetie 
20 geom£trique dans leurs imagettes. 

On peut faire l'hypothese que, pour la 
plupart des images rencontrees en pratique (a 
1' exception notamment des images texturees) , si une 
portion telle que J3 a un contenu insignif iant , 
25 c'est-a-dire qui ne revele rien par rapport aux valeurs 
significatives de l'imagette, les portions homologues k 
plus grande echelle, ici J2 et Jl, seront elles aussi 
insignif iant es . 

Conformement a 1' invention, on commence 
30 done, pour les imagettes des details, a retroprojeter 
les details a 1' echelle la plus petite, puis les 
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details aux echelles de plus en plus grandes. Un seuil 
numerique est applique aux coefficients de l'ondelette, 
c'est-a-dire aux valeurs que prend la transformee dans 
l'imagette consideree. Une valeur inferieure a ce .seuil 
5 donne une portion insignif iante, telle que J3 . 
Cependant, les portions insignif iantes des imagettes ne 
sont pas reconstruites, c'est-a-dire que les calculs de 
retroprojection ne sont pas effectues pour elles. 

On procdde concretement en construisant un 
10 masque numerique avant de retroprojeter l'imagette. Sur 
les details horizontaux Dh, le masque est construit 
pour la premiere fois pour l'imagette Dh3 . II prend une 
valeur egale a 0 pour les portions insignif iantes 
telles que J3 et egale a 1 ailleurs. Les coefficients 
15 du masque se suivent dans un ordre determine, par 
exemple ligne apres ligne. Les calculs de 
retroprojection sont appliques sur l'imagette 
consideree dans 1' ordre des coefficients du masque. 
Quand un coefficient est egal a. 0, aucun calcul n'est 
20 effectue pour le point correspondant de l'imagette, 
auquel on attribue une valeur nulle dans l'imagette 

retropro j ect£e . 

En abordant 1' image suivante des details 
horizontaux (Dh2), le masque numerique permet de ne pas 

25 considerer la portion J2 , dont les points ne sont pas 
trait6s par 1' unite de calcul qui effectue la 
retroprojection. Si par exemple le masque de l'imagette 
Dh3 possede un coefficient nul a la ligne i et la 
colonne j, on prevoit que les quatre points des lignes 

30 2i et 2i+l, et des colonnes 2j et 2j + l de l'imagette 
Dh2 auront aussi des valeurs insignif iantes . Les 
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calculs de r£troprojection ne seront pas effectu£s pour 
ces points, 

Un masque numerique est ainsi construit 
pour chaque niveau de decomposition. Le masque 
5 decrivant l'imagette Dh2 comprendra par convention des 
coefficients 6gaux a 0 pour toute portion homologue a 
une portion aux coefficients nuls (comme J3) du masque 
de l'imagette correspondante & plus petite 6chelle ; 
pour determiner les valeurs 1 ou 0 des autres points du 

10 masque de Dh2, on utilisera encore les comparaisons 
avec le seuil convent i onne 1 . D' autres portions 
insignif iantes, ayant des valeurs nulles de masque, 
pourront apparaitre. C'est ce qu'on a illustr£ sur les 
imagettes Dvl, Dv2 et Dv3, des decompositions 

15 successives des details verticaux. L'imagette Dv3 
comprend une portion insignif iante K3 dont on retrouve 
les homologues K2 et Kl sur les decompositions a plus 
grande gchelle. L'imagette Dv3 ne comprend, dans cet 
exemple, aucune autre portion insignif iante, mais il a 

20 et£ possible d'en trouver trois autres, notees L2 , M2 
et N2, sur l'imagette suivante Dv2 . Au niveau de 
decomposition a plus grande gchelle, celui de 
l'imagette Dvl, on ne s'occupera pas de retroproj ecter 
les points situ£s aux portions LI, Ml et Nl homologues 

25 de L2, M2 et N2 . 

On va maintenant examiner d' autres 
particularites du proced§ de 1 1 invention favorables, a 
une acceleration des calculs. 

La premiere repose sur l'£galite entre la 
30 transform^ de Fourier d*une projection de 1 1 image k un 
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angle fixe et la transformee de Fourier de 1' image sur 
vine droite de meme angle passant par l'origine. 

En se reportant a la figure 3, ou un jeu de 
projections a ete convert i dans le domaine de Fourier 
5 pour donner des projections de nature f requentielle 
dans le systeme d'axes notes ^1 et £2, le support du 
jeu de projection dans le domaine de Fourier comprend 
les valeurs comprises entre -v 0 et v 0 pour £i comme pour 
£ 2 . Une decomposition en ondelette, comme celle de la 
10 figure 1, fait apparaitre 1' approximation encore notee 
AA dans les frequences les plus basses, autour de 
l'origine £l - C,2 = 0, alors que les details se 
trouvent de part et d' autre de cette approximation, a 
des frequences de plus en plus grandes pour les details 

15 a grande echelle. 

On a represent^ une droite passant par 

l'origine, en faisant un angle 6 avec l'axe horizontal 

£l. Cette droite traverse 1' approximation AA, ainsi que 

les imagettes des details verticaux Dvl et Dv2 . Mais on 

20 est assure que les projections faisant cet angle 6 
seront tout a fait inutiles pour les calculs de 
retroprojection des details diagonaux Dd et horizontaux 
Dh puisque la droite d' angle 9 passe a l'ecart de leurs 
imagettes. L' application de 1' invention comprend alors 

25 une selection, pour chacune des categories d' imagettes, 
des angles de projections qui seront utilis6s dans les 
calculs. Le calcul sur le support des projections est 
elementaire . 

La figure 4 reprend la division d'une image 
decomposee en ondelette et transposee dans le domaine 



30 
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de Fourier. Une reconstruction parfaite peut §tre 
obtenue avec l'emploi d'un nombre determine de 
projections, dependant de la fr6quence (discretisation) 
de reconstruction de 1' image de l'objet. Ainsi, pour 
5 : reconstruire 1 ' approximation AA, on peut demontrer 
• qu'il ' suff ira de choisir un nombre de projections, 
parmi celles qui ont ete faites, correspondant a une 
frequence maximale vl servant de rayon a un cercle 
circonscrivant la representation f requentielle de 
10 1' approximation AA dans la decomposition. De meme, les 
details a petite .echelle seront completement rendus en 
utilisant un nombre de projections correspondant a la 
frequence v2 dans le cercle ayant ce rayon, dans le 
systeme d'axes £l, Z>2, et circonscrivant les 
15 representations f requentielle s du groupe de details 
Dv2, Dd2 et Dh2 a meme echelle. Le support des 
projections dans le domaine de Fourier permet done de 
connaitre aisement les frequences maximales necessaires 
pour la r6troprojection parfaite des imagettes 
20 respectives. 

La mise en oeuvre de 1' invention pour 
exploiter ces deux particularity comprendra done, a la 
retroprojection de chacune des imagettes, une selection 
des projections utiles pour cette retroprojection, en 
25 ecartant les autres ; et event uellement une restriction 
du nombre de projections utiles ef fectivement 
exploitees par le calcul pour n'en garder que le nombre 
utile . 

La figure 5 est un organigramme de 
30 1' ensemble du procede decrit ici. 
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REVENDI CATIONS 

1. Procede de reconstruction d'une image a partir de 
jeux de projections de cette image, comprenant 

5 success ivement : 

- une serie de decompositions successives 
des jeux de projections en ondelette dormant des 
imagettes des jeux des projections, comprenant des 
images d'une approximation (AA) et des series 

10 successives (Dd, Dh, Dv) de details homologues de 
chaque jeu, 

- dans chacune des series et success ivement 
pour des imagettes de details ayant un nombre croissant 
de points, une recherche de portions insignif iantes 

15 jugees pour depourvues de contenu, et une recherche de 
portions homothetiques des portions insignif iantes dans 
chacune des imagettes de details qui suivent dans la 
serie ; 

- des ret roproject ions des imagettes des 
20 jeux de protection en imagettes de 1 ' image a 

reconstruire, avec une omission des retroprojections 
pour toutes les portions insignif iantes et toutes les 
portions homothetiques, 

et une combinaison des imagettes a 
25 reconstruire par une inversion de decomposition donnant 
ladite image. 

2. Procede de reconstruction d'une image suivant la 
revendication 1, caracteris6 en ce qu'il comprend une 
selection de regions d' angles (9) des jeux de 
30 projections qui sont utilises dans les 
retroprojections . 
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3. Procede de reconstruction d'une image suivant la 
revendication 2, caracterise en ce qu'il comprend une 
selection d'un nombre de projections qui sont utilisees 
dans la rStroproj ections . 

4. Proced6 de reconstruction d'une image suivant la 
revendication 2, caracterise en ce que la selection des 
regions d' angle est real i see en fonction d'un support 
d'une transformee de Fourier des jeux de projections 
decomposes en ondelette. 

5. Procede de reconstruction d'une image suivant la 
revendication 3, caracterise en ce que la selection du 
nombre de projections est realis6e en fonction de 
frequences maximales d'un support d'une transform6e de 
Fourier des jeux de projection decomposes en ondelette. 
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MESURES (JEUX BIBIMENSIONNELS 
DE PROJECTIONS, OU PREMIERE 
IMAGE DE L'OBJET) 



REAMENAGEMENT DES 
MESURES CONIQUES EN 
MESURES EN EVENTAILS- 
PARALLELES OU DES 
MESURES EN EVENTAIL 
EN MESURES EN PARALLELES 



DECOMPOSITION EN ONDELETTE 
DE L' IMAGE EN IMAGETTES 



SELECTION DES 
PROJECTIONS UTILES 
POUR CHAQUE IMAGETTE 



IMAGETTES DES DETAILS A 
LA PLUS PETITE ECHELLE 



DETERMINATION D'UN 
MASQUE DES IMAGETTES 



FIG. 5 



RETROPROJECTION 
DES IMAGETTES DE 
MEME ECHELLE. 



PASSAGE A 
L'ECHELLE 
SUPERIEURE 




NON 



SYNTHESE DE L'IMAGE 
DE L'OBJET PAR COMBINAISON 
DES IMAGES RETROPROJETEES 
EN ONDELETTE INVERSE 
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